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We present a simple device architecture that allows all-electrical detection of plasmons in a 
graphene waveguide. The key principle of our electrical plasmon detection scheme is the non¬ 
linear nature of the hydrodynamic equations of motion that describe transport in graphene at room 
temperature and in a wide range of carrier densities. These non-linearities yield a dc voltage in 
response to the oscillating field of a propagating plasmon. For illustrative purposes, we calculate 
the dc voltage arising from the propagation of the lowest-energy modes in a fully analytical fashion. 

Our device architecture for all-electrical plasmon detection paves the way for the integration of 
graphene plasmonic waveguides in electronic circuits. 


Introduction .—The two-dimensional (2D) electron liq¬ 
uid in a doped graphene sheet ^ supports plasmons with 
energies from the far-infrared to the visible, depending 
on carrier concentration^. Although they share similar¬ 
ities with plasmons in ordinary parabolic-band 2D elec¬ 
tron liquids^, plasmons in graphene are profoundly dif¬ 
ferent. From a fundamental point of view, their disper¬ 
sion relation is sensitive to many-body effects even in the 
long-wavelength limit^. More practically, plasmons in 
graphene are easily accessible to surface-science probes 
and opto-electrical manipulation since they are exposed 
and not buried in a quantum well. 

Plasmons in graphene are also substantially different 
from those in noble metals. Indeed, recent near-held opti¬ 
cal spectroscopy experiments®^® have demonstrated that 
plasmons in graphene display gate tunability and ultra¬ 
strong held conhnement. Moreover, low damping rates 
can be achieved by employing graphene samples encapsu¬ 
lated in hexagonal boron nitride thin slabs®^^®. For these 
reasons, graphene plasmonics has recently attracted a 
great deal of interest^®. Graphene plasmons may al¬ 
low for new classes of devices for single-plasmon non- 
linearities^^, extraordinarily strong light-matter inter¬ 
actions^®, deep sub-wavelength metamaterials^®^^®, and 
photodetectors with enhanced sensitivity®®’^^. 

A key ingredient of a disruptive plasmonic platform is 
the ability to efhciently detect plasmons in all-electrical 
manners. Some progress has been made in this direction 
in conventional noble-metal-based plasmonics. Falk et 
for example, were able to couple plasmons in Ag 
nanowires to nanowire Ge held-effect transistors. Built-in 
electric helds in the latter are used to separate electrons 
and holes before recombination, thereby giving rise to a 
measurable source-drain current. Similarly, Neutens et 
al.^^ employed an integrated metal-semiconductor-metal 
detector in a metal-insulator-metal plasmon waveguide. 

While graphene plasmons have been detected and 
studied in a multitude of ways^®, including electron 
energy loss spectroscopy®^, polarized Fourier transform 



FIG. 1. (Color online) Schematics of our electrical plasmon 
detector. A graphene strip of width W is encapsulated be¬ 
tween two dielectrics (semi-transparent slab above and dark 
green slab underneath graphene). A back gate (dark blue 
slab), separated by a distance d from the graphene sheet and 
held at a voltage Vg, is used to control the average carrier 
density h in graphene. At one end of the strip, a plasmon is 
launched by using e.g a metallized atomic force microscope 
tip illuminated by lightDue to non-linearities in the hy¬ 
drodynamic equations, a dc electrical potential difference 5V 
is measured between probe electrodes placed at positions ri, 
r 2 , and rs and a reference electrode placed at the other end 
of the strip. The quantity 5V provides a direct measurement 
of the ac electric field of a propagating plasmon. 


infrared spectroscopy^®, and near-held optical spec- 
troscopy®^®, a protocol for all-electrical detection of these 
modes is still lacking. 

In this work we present a device architecture that 
allows all-electrical detection of plasmons in graphene 
waveguides. In our scheme, all-electrical detection is not 
enabled by the integration of a detector in a graphene 
plasmon waveguide (GPW) but rather by the intrinsic 
non-linear terms in the hydrodynamic equations that de¬ 
scribe transport in the 2D massless Dirac fermion (MDF) 
liquid^ hosted by graphene. Non-linearities enable the 
emergence of a rectihed (i.e. dc) component 6V{r) of 
the ac electric held of a propagating plasmon, which can 
be measured by a suitable geometry of ohmic contacts 
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placed along the GPW, as shown in Fig. 1. We now 
present a calculation of the spatially-dependent electri¬ 
cal signal SV{r). 

Hydrodynamic theory. —We consider a GPW with 
transverse (longitudinal) size W {L with L ^ W), which 
is embedded between two insulators with dielectric con¬ 
stants ei (above the GPW) and 62 (below the GPW). 
Here, “longitudinal” and “transverse” refer to the plas- 
mon propagation direction— x in Fig. 1. 

We would like to describe ac transport in a GPW by 
employing the theory of hydrodynamics^^. We there¬ 
fore need to assess whether experimentally relevant re¬ 
gions of parameter space exist in which this theory is 
applicable. First, at room temperature and for typical 
carrier densities (h ~ 10^^ cm“^-5 x 10^^ cm“^), the 
mean-free-path £ee = vfTbb for electron-electron colli¬ 
sions in graphene is short^®’^^, i.e. £ee — 100-150 nm. 
Here, vf ~ 10® m/s is the graphene Fermi velocity^® and 
Tee — 100 fs = 10 “^® s is the electron-electron collision 
time^®’^^. Second, for hydrodynamics to provide a cor¬ 
rect description of the response of the system at finite 
frequencies, it must also be uiTee 'C 1 , where ui is the 
external-excitation angular frequency. The value of Tee 
given above constraints the maximum external-excitation 
frequency to be /max = l/(27rTee) < 3 THz. We there¬ 
fore conclude that, for n ~ 10^^ cm“^-5 x 10^^ cm“^, 
w < 27r/max, and T = 300 K, transport in GPWs with 
characteristic dimensions L,W £ee is accurately de¬ 
scribed by hydrodynamic equations of motion^®. Related 
continuum-model descriptions of plasmons in GPWs have 
been employed in Refs. 29-32. 

The set of hydrodynamic equations consists of i) the 
continuity equation, 

dtn{r,t)■[n{r,t)v{r,t)\ = Q , (1) 

and ii) the Navier-Stokes equation^® 

mcn{r,t)Dtv{r,t) = —en{r,t)E{r,t)+r]V^v{r,t) . (2) 

In Eqs. (l)-(2), n{r,t) is the carrier density and v{r,t) 
is the drift velocity. In Eq. (2), nic = h\^TT^fvF is the 
graphene cyclotron mass^®, with n = CVaje the aver¬ 
age electron density and Vq the back-gate voltage (see 
Fig. 1), and Dt = dt + v{r, t) • V is the convective deriva¬ 
tive^®. The electric field E{r,t) = —V^(r,t) is the gra¬ 
dient of the electrostatic potential ^{r,t) (we neglect re¬ 
tardation effects). Finally, rj is the shear viscosity of the 
2D electron liquid®’^®. For future purposes, we also in¬ 
troduce the kinematic viscosity^® 


It can be shown®® that, in the hydrodynamic ujTee 1 
limit, 1 / ~ WpTee/4. With the values of vf and Tee given 
above, we find ly ~ 250 cm^/s. In writing Eq. (2) we 
have neglected a term due to the bulk viscosity C since 
this quantity vanishes at long wavelengths®’^®. 


We highlight two non-linear terms in Eqs. (l)-(2): 
a) the non-linear coupling between n{r,t) and v(r,t), 
which is present in Eq. (1), and b) the non-linear term 
\v{r,t) ■ W]v{r,t) in Eq. (2), representing the convective 
acceleration^®. 

Momentum-non-conserving collisions, such as those 
due to the friction of the electron liquid against the disor¬ 
der potential, can be taken into account phenomenolog¬ 
ically by adding a term of the type —mcjn{r,t)v{r,t) 
on the right-hand side of Eq. (2), where 7 is a damp¬ 
ing rate®^. Furthermore, corrections to Eq. (2), stem¬ 
ming from the pseudo-relativistic nature of MDF flow in 
graphene, can be easily incorporated into the theory®®’®® 
and have been demonstrated to yield stronger rectified 
signals®®. 

Finally, to close the set of equations, we need a rela¬ 
tion between d)(r,t) and n{r,t). This depends on the 
screening exerted by dielectrics and conductors near the 
GPW. If a metal gate is positioned underneath the GPW 
at a distance d <C W, k~^, where k is the plasmon wave 
vector, the following local relation exists®®: 

d>(r,t) ft!-^dn(r,t) , (4) 

where C = e 2 /( 47 rd) is a capacitance per unit area and 
6 n{r, t) = n{r, t) — n. Eq. (4) greatly simplifies the theo¬ 
retical analysis and, in fact, allows us to solve the prob¬ 
lem in a fully analytical fashion®^, as we now detail. 

Eqs. (l)-(4) need to be accompanied by boundary con¬ 
ditions. As explained in Appendix A and in Ref. 38, 
we fix Vy{x,y = 0,VF) = 0 and dxVy{x,y = 0,W) + 
dyVx{x,y = 0,W) = 0. 

Linear response theory and plasmons .— The GPW sup¬ 
ports collective charge density oscillations, i.e. plas¬ 
mons®, which propagate along the x direction and are 
confined in the y direction. To calculate the frequency 
spectrum and potential profiles of these modes we have 
to linearize Eqs. (I)-(2) and (4). We write n{r,t) = n + 
ni{r,t) + n2{r,t) + ..., v{r,t) = Vi{r,t)+V2{r,t) + ..., 
and 4)(r,t) = <i)i(r,t) -|- 4 ) 2 (r, t).... Here ni{r,t), 
vi{r,t), and ^i{r,t) [n2{r,t), V2{r,t), and <i) 2 (r,t)] de¬ 
note first-order [second-order] corrections with respect to 
equilibrium values (by “equilibrium” we here mean the 
state of the GPW in which a plasmon is not propagat¬ 
ing). In the linearized theory we retain only terms of 
the first order. All the relevant details are reported in 
Appendix B and G. 

For the sake of simplicity, we assume a uniform equi¬ 
librium electron density in the GPW, disregarding the 
well known inhomogeneous doping n —>■ n(y) that arises 
due to a back gate. Plasmons in back-gated waveguides, 
however, have been demonstrated®® to be similar to those 
of uniformly doped waveguides, provided that the Fermi 
energy is appropriately scaled to compensate for the sin¬ 
gular behavior of the carrier density n{y) as y —>■ 0, W. 

Plasmon modes are labelled by a wave number k (stem¬ 
ming from translational invariance along the x direction) 
and a discrete index n = 0,1,2,.... The associated ac 
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electrical potential is given by 

where 

1 J 1, for n = 0 
^/W I •\/ 2 cos[n 7 r 2 //t 4 ^], for n 0 

The mode dispersion reads as following 


(5) 


( 6 ) 




where = irnlW, and s = ^/e‘^hf{Cmc) 

is the hydrodynamic speed of sound. It is useful to intro¬ 
duce the following natural frequency scale: flo = sjW. 
Setting external ( 7 ) and internal (v) dissipation to zero 
in Eq. (7) we find the expected result a;„(fc) = sKn- The 
lowest-energy n = 0 mode shows an acoustic dispersion 
due to screening by the back gate. Modes with n ^ 0 are 
gapped, i.e. a;„(fc —>■ 0) = mrflo- The fundamental fre¬ 
quency is a;„=i/(27r) = flo/2 — 1.0 THz for IT = 3 ^m, 
d = 100 nm, 62 = 3.9, and fi = 10^^ cm“^. Dispersion 
relations and mode profiles for the above set of param¬ 
eters are shown in Fig. 2. In the approximation (4) the 
results do not depend on ei. 

When the n-th eigenmode of the GPW is excited by 
an external perturbation with frequency w, it propagates 
with a complex wave number 


fcn(w) = 


-I- lUJ-J „ 

W-^- 

— lULU 


( 8 ) 


The wave number, 3?e(fc„), and inverse damping ratio, 
7 “^ = ^e{kn)/‘^m{kn), of the launched plasmon depend 
only on the excitation frequency w and not on details of 
the tip-sample coupling®^®. Physically, the dimension¬ 
less number 7 “^ controls the plasmon extinction length 
(.n = l/3m(fc„) = 7“U„/(27r), with A„ = 27r/3fie(fc„) 
the plasmon wavelength. With the value of v given above 
and 7 = 0.3 flo, the inverse damping ratio of the n = 0 
mode is ~ 23, while 7 f ^ ~ 7 for n = 1. 

Second-order theory of all-electrical detection .—The 
rectified signal can be calculated by keeping track of 
the second-order terms n 2 {r,t), V 2 {r,t), and ^ 2 {r,t) in 
the expansion of the hydrodynamic variables. Physically, 
the second-order response describes interactions between 
propagating modes. If only one mode propagates, a dc 
signal due to self-mixing of the plasmon field is gener¬ 
ated. If more than one mode propagates, also interference 
terms will be generated. On general grounds^^, we expect 
that the second-order response is composed by an oscil¬ 
lating component at frequency 2uj (i.e. second-harmonic 
generation) and a steady component. Since we are inter¬ 
ested in detecting a dc signal, we can extract the rectified 
voltage 8 V (r) from the time average over one period of 
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FIG. 2. (Color online) Panel (a): Dispersion relation /n(fc) = 
ujnik)/(fin) of four low-energy plasmon modes (n = 0,... ,3) 
in a GPW with the following parameters: IT = 3 /rm, 
d = 100 nm, £2 = 3.9, n = 10^^ cm“^, 7 = 0.3 Do, and 
n = 250 cm^/s. Solid (dashed) lines correspond to even 
(odd) modes. Panels (b)-(e): Corresponding electrical po¬ 
tential profiles <l?i(r,t) evaluated at x = 0 and t = 0 and 
plotted as functions of the transverse y G [0, IT] coordinate. 
These results have been obtained by imposing electron density 
fluctuations equal to"*® 3n/n = 1%. 


the external radiation of the second-order potential fluc¬ 
tuations: 6V{r) = ( 4 ) 2 (r,t)). Averaging over time the 
second-order equations as explained in Appendix B, we 
obtain 

■ Sv{r) =■ {ni{r,t)vi{r,t)) (9) 

n 

and 

— —VSV{r) -\- 'ySv{r) — vV'^5v{r) 

me 

= -e{[vi{r,t) ■ \/]vi{r,t)) - ^(ni(r,t)V^'Ui(r,t)) , 

n 

( 10 ) 

where 5v = {v 2 {r,t)) is the time average of the veloc¬ 
ity fluctuations. We urge the reader to note that in 
Eq. (10) we have introduced a dimensionless parameter, 
e, which allows us to keep track of the role of differ¬ 
ent non-linearities in determining the rectified signal. By 
setting e = 0 one neglects the convective non-linearity in 
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the Navier-Stokes equation. Moreover, by setting e = 0 
and ^ = 0, the Navier-Stokes equation reduces to the 
linearized Euler equation^^, which leads to the standard 
Drude formula for the local conductivity. However, a fi¬ 
nite rectified signal 5V{r) exists in this case too and is 
entirely due to the non-linear n(r,t)v(r,t) coupling in 
the continuity equation. 

Eqs. (9)-(10) are crucial since they relate the second- 
order quantities SV and Sv to the quantities ni{r,t) and 
Vi{r, t), which have been calculated in the linearized the¬ 
ory. Furthermore, they can be used to calculate the dc 
signal 5V(r) in response to plasmon propagation in any 
desired geometry. As stated above, SV{r) can be mea¬ 
sured by employing a set of ohmic contacts as in Fig. 1. 

We now evaluate SV{r) for the experimentally rele¬ 
vant case in which plasmons are launched at a specific 
location r* = {0,y*) with y* e [0,W] in the GPW. The 
quantity 5V{r) for x ^ W and arbitrary y can then be 
calculated according to the following procedure, which is 
typical of a scattering problem, a) For x^W the plas¬ 
mon velocity field can be written as a sum over prop¬ 
agating modes (i.e modes with a;„(fc) < w, where uj is 
the angular frequency of the stimulus that launches plas¬ 
mons). All the other modes, which can be excited near 
r*, exponentially damp out at large distances since they 
have a purely imaginary k —see Eq. ( 8 ). Furthermore, 
as shown in Appendix C, the plasmon velocity field Vi is 
irrotational at large distances, i.e. V x = 0 for x W. 
Since V x V(j){r) = 0 for an arbitary scalar function 4){r), 
we conclude that Vi{r,t) for a; ;g> W can be written as 
the gradient of a scalar function. In the language of scat¬ 
tering theory, we have built the so-called asymptotic so¬ 
lution, which we denote by v[^\r,t). b) Let us imagine 
that an external perturbation with frequency uj launches, 
for example, an arbitrary linear combination with com¬ 
plex coefficients of the n = 0 and n = 1 GPW modes. 
Because of a), we can write the corresponding asymptotic 
velocity field as 

v[^\r,t) = ^V[(l - 

-k -k c.c. , ( 11 ) 

where a: > 0 and the functions <Pniy) have been intro¬ 
duced earlier in Eq. ( 6 ). In Eq. (11) A = is an 

unknown amplitude (here v has physical dimensions of 
a velocity), which can be estimated as discussed below, 
kn = 3?e[A:„(cLi)] and /3„ = 'ism[kn{uj)] with fc„(w) as in 
Eq- ( 8 ), ^ G [ 0 , 1 ] is a real parameter that allows us to in¬ 
terpolate between the case in which only the n = 0 mode 
is launched (^ = 0 ) and the case in which only the n = 1 
mode is launched (^ = 1 ), and e*“ (with a real) is the rel¬ 
ative phase between the two modes. In the case in which 
plasmons are launched at r* by using a metallized tip 
illuminated by light®^®, ^ depends on the tip-sample cou¬ 
pling: for example, for a tip placed symmetrically with 
respect to the GPW axis ^ vanishes. In practice, ^ can be 
varied by moving the tip along the y direction, c) With 
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FIG. 3. (Color online) Panel (a) The dc potential 5V{r) (in 
millivolts) as a function of x (in /rm), calculated at the edges 
of the GPW, i.e. at y = 0 (solid, dashed, and dotted lines) 
and y = W (dash-dotted line). These numerical results have 
been calculated by setting the following parameter values in 
Eq. (11): ^ = 0 and a = 0 (solid line), ^ = 1 and a = 0 
(dashed line), ^ = 1/2 and q = 7 r/4 (dotted and dash-dotted 
lines). The other parameters are as in Fig. 2. Finally, we have 
taken tj/(27r) ~ 1.2 THz in Eq. (8). Note that ijj/{2'k) is 20% 
larger than the fundamental frequency for the laser to be able 
to excite the two lowest modes of the GPW. The scale is in 
millivolts. In all cases the value of A has been chosen to yield 
5n/n = 1%. This normalization implies different values of A 
for different values of a. Panel (b) Spatial map of the dc 
potential 5V{r) calculated by setting f = 1/2 and a = 7 r/4. 
The color bar shows the potential amplitude in millivolts. All 
other parameters are as in panel (a). 


the velocity field in Eq. (11), one can easily calculate the 
asymptotic density profile n^i\r,t) from the continuity 
equation, d) The quantities n^i\r,t) and v^\r,t) are 
then used to calculate the temporal averages that appear 
on the right-hand side of Eqs. (9)-(10). e) Finally, SV{r) 
is found by solving Eqs. (9)-(10). 

Simple and compact analytical expressions, obtained 
by following steps a)-e), are available for SV{r) in the 
extreme cases C = 0 and ^ = 1 and are presented Ap¬ 
pendix D. In the general ^ ^0,1 case an oscillatory term 
with spatial periodicity 27r/(fco — fci) appears along the x 
direction due to interference of the two modes in Eq. (11). 
Illustrative numerical results can be found in Fig. 3. No¬ 
tice that the dc signal is < 0.5 mV and that its spatial 
extension is ~ 20 /rm. The dc potentials on the top and 
bottom GPW edges are not equal in the case ^ yf 0,1 
since Eq. (11) is a superposition of modes with different 
parity. The quantity A in Eq. (11) was estimated with 
reference to Refs. 5-9, where a metallized tip is used to 
launch plasmons. In this case, A can be calculated start- 
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ing from the amplitude of density oscillations Sn in units 
of n, created by the tip at r*. The results in Fig. 3 have 
been produced by using"^'^ Sn/n = 1% at a; = 0. This 
“normalization” condition yields different values of A for 
different values of a. In other words, A measures how 
well the tip couples to the linear combination of modes in 
Eq. (11). The ^ = 0 mode has a better coupling to the tip 
(and therefore yields a larger dc signal) than the linear 
superposition of n = 0 and n = 1 modes corresponding 
to ^ = 1/2 and a = 7 r/ 4 . 

Conclusions .—In summary, we have discussed an ar¬ 
chitecture based on a graphene waveguide where elec¬ 
trical detection of plasmons may be experimentally 
achieved. We have theoretically demonstrated that rec¬ 
tification of the ac field of a propagating plasmon, which 
is enabled by non-linear terms in Eqs. (l)-(4), yields a 
spatially-dependent dc signal SV{r). The experimental 
exploitation of similar non-linearities has recently led^'^ 
to room-temperature graphene THz photodetectors. We 
stress that SV{r) can be calculated from Eqs. (9)-(10) 
and can be measured by using lateral probe contacts as 
in Fig. 1. Simple analytical expressions for 6 V{r) have 
been given in Appendix D for the cases ^ = 0,1. Nu¬ 
merical results for the general case ^ yf 0,1 have been 
presented in Fig. 3. Such values of dc voltages can be 
easily measured. 
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Appendix A: Boundary conditions 

Hydrodynamic equations need to be accompanied by 
appropriate boundary conditions (BCs) at the bound¬ 
aries of the region occupied by the electron fluid. In the 
general case of viscous flow, we need BCs on the normal 
and tangential components of the current evaluated at 
the boundaries of the GPW. 

The BC on the normal component can be simply de¬ 
rived by requiring that no charge exit the boundary of the 
system. The normal component of the current J = nv 
evaluated at the boundaries must therefore vanish 

hiJi = nhiVi = 0 , (Al) 

where n is a unit vector normal to the boundary. In our 


simple rectangular geometry this translates to 

'^yi^Ty)\y=0^W ~ ® ■ (-^ 2 ) 

In the description of the motion of liquids near a fixed 
surface, the following BC, which physically corresponds 
to the vanishing of the tangential component of the ve¬ 
locity field, is usually assumed^^: 

eijfiiVj = 0 . (A3) 

where is the 2D fully-antisymmetric tensor. Eq. (A3) 
can be justified for molecular liquids^^, where interac¬ 
tions between two molecules of the liquid are essentially 
of the same type of the interaction between a molecule 
in the liquid and one in the surface of the container. 

For the electron liquid in graphene, however, in which 
boundaries are defined by etching or electrostatic gat¬ 
ing, we cannot identify a mechanism through which the 
boundary can exert a tangential force on the liquid. For 
this reason we chose to apply the following BC^®: 

= 0 ■ (A4) 

Physically, the previous equation represents the van¬ 
ishing of the tangential force on the liquid, and is 
commonly used^® for liquids near a free surface. In 
Eq. (A4) a'ji, is the 2D viscous stress tensor, i.e. cr'^, = 
V (djVk + dkVj — Sjkdivi). In writing the previous expres¬ 
sion for a'jf. we have set to zero a contribution due to the 
bulk viscosity 

In our geometry, Eq. (A4) reduces to: 

[d^Vy{x, y) + dyvffx, = 0 , (A5) 

as reported in the main text. 

Appendix B: Perturbative solution of the 
hydrodynamic equations 

In the main text we have solved to the coupled non¬ 
linear equations (l)-(2) by employing the following per¬ 
turbative ansatz^^: 

n{r,t) = n + Xni{r,t) + X^n 2 {r,t) -|- ... 

v{r,t) = Xvi{r,t) + X^V 2 {r,t) + ... , (Bl) 

^{r,t) = A<i)i(r,t) -I- X'^^ 2 {r,t) -I- ... 

where A is a bookkeeping parameter, which will be set to 
unity at the end of calculations. 

Inserting the above expressions into Eqs. (l)-(2) of the 
main text and retaining only terms of the first-order in 
A, we obtain the linearized continuity equation, 

dtni{r,t)+X7 ■ [nvi{r,t)] = 0 , (B2) 

and the linearized Navier-Stokes equation, 

dtVi{r,t) = —V4>i(r,t) - jvffrff) + . 

me 

(B3) 
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These two equations have been used in the main text to calculate the spectrum of plasmons in the GPW. Note that, 
if the last equation is written in Fourier transform with respect to both space and time, it leads to a conductivity of 
the form 


cr{q,U}) 


i'DolTT 

w + Z7 + ’ 


(B4) 


where I?o = irne'^/mc is the (non-interacting) Drude weight^. Viscosity yields a non-local correction to the ordinary 
Drude formula. 

At second order in A we obtain the following equations: 


dtn 2 {r,t) + V • [nv 2 {r,t) + ni{r,t)vi{r,t)] = 0 


(B5) 


and 


rricn [dtV2{r-t) + {vi{r;t) ■ V)ui(r;t)] -|- mcni{r;t)dtVi{r;t) = enV 4 > 2 (r;t) -|- eni(r; t)V 4 >i(r; t) 
-mcjhV2{r-, t) - TOc7ni(r; t)vi{r-, t) + ri'V'^V2{r; t) . 

The latter can be simplified thanks to the linearized Navier-Stokes equation (B3), leading to: 

dtV2{r,t) -^V 4 > 2 (r,t) +'YV2{r,t) - vV‘^V2{r,t) = V^t>i(r, t) - e[ui(r,t) • V]ui(r,t) . 

rric n 

Averaging over time both sides of Eqs. (B5) and (B7) we finally find Eqs. (9)-(10) of the main text. 


(B6) 


(B7) 


Appendix C: Normal modes of oscillation: plasmons 

In this Section we use the linearized hydrodynamic equations (B2)-(B3) to show that the GPW supports collective 
charge density oscillations, i.e. plasmons. 

We look for oscillating solutions where all the quantities vary in time as f{r,t) = Re(*’)]) where uj is, in 
general, a complex number (to allow for modes with a finite linewidth). Eqs. (B2)-(B3) then transform into: 

- + V ■ [hvi^^{r)] = 0 (Gl) 


and 

e2 

- iujvi^^{r) = -7:;Vni,;^(r) - , (G2) 

rricC 

where the electric potential <i>i has been rewritten by emptying the local capacitance approximation, i.e. Eq. (4) of 
the main text. 

We can eliminate in Eq. (G2) by using Eq. (Gl). We find 

e2 

w^Vi,c;(r)H-=0. (G3) 

rricC 

We can take advantage of translational invariance along the x direction by introducing the following ansatz for 


ViAr) = 


oo oo 

x'^an(Pn{y) +y'^bn1pn{y) , 

n—0 n—1 


(G4) 


where the functions (pn have been introduced in the main text, while ijjn{y) = \/2/IT sin(n7ry/IT). It can be easily 
checked that {(p„(y),n = 0,1,... } and {^„(y),n = 1,... } represent complete sets of orthonormal functions in the 
interval y S [0, IF]. Furthermore, the ansatz (G4) respects the BGs (A2) and (A5). 

Substituting Eq. (G4) in Eq. (G3) and using standard completeness theorems, we can transform Eq. (G3) into a 
set of decoupled matrix equations, one for each n: 


-I- iyw -I- + n^7r^/IF^)] 


hn 


^ f —ikWn-K 

^ ikWmr 



(G5) 
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where fig = e^n/{CrricW^) is the square of a natural frequency scale, as in the main text. 

The 2x2 matrix on the right-hand side of Eq. (C5) has the following eigenvectors 

The solutions are then given by Eq. (C4) where, for every n’ = 0,1,..., the vector (a„, 6„) = 5n,n'^h O'" (“n, bn) = 

^n,n' ^ n’ 

The solutions built by using €>„ in Eq. (C6), which has eigenvalue k'^W^ + lead to an irrotational spatial 

profile of the type: 


vi^a,{r) = [ikipn{y)x - n'K'iljn{y)v\= AV [e^’^'^ifniy)] , (C7) 

with uj = LOn{k) given by the dispersion relation in Eq. (7) of the main text. 

The second eigenvector, fl/n, with zero eigenvalue, leads to solenoidal solutions with no density fluctuations. These 
modes have a purely imaginary frequency (because they do not experience a restoring force that sustains them) 
and cannot therefore propagate along the GPW. Even if they are excited by an external perturbation (note that an 
external scalar potential does not couple to these modes) they decay exponentially for a: ^ VE. We can therefore 
safely exclude them from the construction of the asymptotic solution v^^\r^t) in the main text. 


Appendix D: Analytical results for the rectified signal 

We here report our main analytical results for the rectified signal in the extreme cases ^ = 0 and ^ = 1—see Eq. (11) 
in the main text. For ^ = 0 we find 


SV(r) = (feg + Pi) + fco7 + (fcg - 3/3g)] 

^ ^ 4/3ow 


(Dl) 


while, for ^ = 1, 

s:^r( _ V I + /3ifci7) (/3? + fc? + gi) + ^Piki (fcf - 3/3f + qt- 2klPl - QqlPl + 2qlkl) 

0V{r) i/e tv j 

¥> 2 ( 2 /) [ew(/3i - ql) + Piki-f] {Pl + kl- ql) + uPiki (fcf - 3/3f - 3gf - 2klj3l + QqlPf + 6qfkj) \ 

72 Ml3!-ql) j ■ 

(D2) 


Here V = je = mcv"^(e has physical di¬ 

mensions of a voltage. Once again, we hasten to stress 
that a finite value of 5V{r) exists also in the ordinary 
Drude iz = e = 0 limit, as it can be explicitly checked 


from Eqs. (D1)-(D2). In this case the rectified signal is 
due to the hydrodynamic n{r,t)v(r,t) non-linear cou¬ 
pling in the continuity equation. 
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